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ABSTRACT 

We show that the period clustering of anomalous X-ray pulsars (AXPs) and 
soft gamma-ray repeaters (SGRs), their X-ray luminosities, ages and statistics 
can be explained with fallback disks with large initial specific angular momentum. 
The disk evolution models are developed by comparison to self-similar analytical 
models. The initial disk mass and angular momentum set the viscous timescale. 
An efficient torque, with (1 — uJl) dependence on the fastness parameter cu* leads 
to period clustering in the observed AXP-SGR period range under a wide range of 
initial conditions. The timescale t for the early evolution of the fallback disk, and 
the final stages of fallback disk evolution, when the disk becomes passive, are the 
crucial determinants of the evolution. The disk becomes passive at temperatures 
around 100 K, which provides a natural cutoff for the X-ray luminosity and defines 
the end of evolution in the observable AXP and SGR phase. This low value for 
the minimum temperature for active disk turbulence indicates that the fallback 
disks are active up to a large radius, ^ 10 12 cm. We find that transient AXPs and 
SGRs are likely to be older than their persistent cousins. A fallback disk with 
mass transfer rates corresponding to the low quiescent X-ray luminosities of the 
transient sources in early evolutionary phases would have a relatively lower initial 
mass, such that the mass-flow rate in the disk is not sufficient for the inner disk to 
penetrate into the light cylinder of the young neutron star, making mass accretion 
onto the neutron star impossible. The transient AXP phase therefore must start 
later. The model results imply that the transient AXP/SGRs, although older, 
are likely to be similar in number to persistent sources. This is because the 
X-ray luminosities of AXPs and SGRs are found to decrease faster at the end 
of their evolution, and the X-ray luminosities of transient AXP and SGRs in 
quiescence lie in the luminosity range of X-ray cutoff phase. Taking the range of 
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quiescent X-ray luminosities of transient AXPs and SGRs ~ 10 33 — 10 34 erg s -1 , 
our simulations imply that the duration of the cutoff phase is comparable to the 
lifetime in the persistent phase for a large range of initial conditions. 

Subject headings: pulsars: individual (AXPs) — stars: neutron — X-rays: bursts - 
accretion, accretion disks 
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INTRODUCTION 



Anomalous X-ray pulsars (AXPs ' 



Woods Thompson 



2006 



Mereghetti 



and soft gamma-ray repeaters (SGRs) (see 



20081 . for a review) are rapidly spinning down young 



neutron stars with persistent X-ray luminosities L x ~ 10 34 — 10 36 erg s -1 in excess of their 
spin-down power. The spin periods of this small population with ~ 15 members is clustered 
between P = 2 — 12 s; placing them in the top right corner of the P — P-diagram, populated 
by only a few radio pulsars. 

The recognition of AXPs and SGRs among the new classes of young objects, along 
with other young neutron star families which do not manifest themselves as radio pulsars, 
like the dim isolated neutron stars (DINs) and the compact central objects (CCOs), 
suggests alternative evolutionary paths for young neutro n stars. The magnet ar model 



( iDuncan fc Thompson 



1992 



Thompson fc Duncan 



19951 ) attributes SGR and AXP bursts 



and other properties to magnetic field stronger than quantum critical field, B c = 4.4 x 10 13 G. 



The presence of high magnetic field radio puls ars wit 



this 



Deing the only criteri on. An alternative (Alpar 



20001 ) or complementary (lEksi &: Alpar 



2003 



i super-critical inferred fields challenges 



2001 



Chatteriee. Hernquist fc Narayan 



Ertan &: Alpar 



20031 ) view suggests the 



presence of remnant disks around young neutron stars left over from the supernova explosion 
as a necessary ingredient causing the diversity. The initial mass and angular momentum 
of the disk present a new set of parameters, in addition to the initial period and dipole 
magnetic moment of the neutron star, that affect evolutionary paths and observational 
classes of neutron stars. 

AXPs and SGRs show bursts with luminosities well above the Eddington limit. This 
cannot be explained by gas accretion onto the neutron star. The magnetar model explains 
these bursts as instabilities in the crust due to strong magnetic stresses and reconfiguration 
of the magnetosphere, while it accounts for the observed spin-down torques of AXP and 
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SGRs by magnetic dipole radiation. In the magnetar model it is hard to explain the period 
clustering of AXPs and SGRs within the f ramework that the per sistent X-ray luminosity is 



a consequence of the magnetic field decay (jPsaltis &: Miller 



20021 ) 



The detection of a debris disk around the AXP 4U 0142+61 (IWang. Chakrabarty fc Kaplan 



20061 ) could be the clue for what makes AXPs and SGRs different from radio pulsars in 



their persistent phases (jErtan. Erkut. Ek§i fc Alpar 



20071 ). The presence of fallback disks 



can address the period clustering of AXPs and SGRs as a result of the long-term interaction 
of the disk with the magnetic dipole field of the neutron star, converging to "equilibrium" 
periods. While the period clustering generated by the earlier models is in agreement with 
observations, there were difficult ies in explaining the spin-periods together with X-ray 



luminosities (jEksi fc Alpar 



2003|) 



In this paper, we have identified three crucial aspects of fallback disk evolution: (1) The 
magnetic torque on a "mildly" fast neutron star in the regime of accretion with spin-down , 



Ertan fc Erkut 



has th e dependence N oc (1 — u>l) on the fastness parameter a;* (see also 
( 120081 )). (2) The disk timescale t Q that determines the initial evolution of the system is not 
arbitrary, but instead set by the initial mass and angular momentum of the disk. (3) The 
outer regions of the disk become passive at temperatures T < T p ~ 100 K. This results in 
the mass-inflow rate M of the disk being turned off rather sharply at ages ^ a few 10 4 — 10 5 
y. In § 2, we present the disk equations and self-similar models for fallback disk evolution. 
The timescale to, which is particularly important for early disk evolution is related to the 
initial disk mass and angular momentum. We also briefly discuss the assumptions and 
results of earlier analytical work on the long-term evolution of fallback disks. Comparison 
of analytical and numerical models is used to set the stage for the numerical calculations 
detailed in § 3. This section also presents the torque model employed in our calculations 
and the justifications for this choice of torque on the basis of work on transient AXPs. 
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The late evolution of the fallback disk is also addressed in 5 3. Fallback disk mode 



in evolutionary calculations so far (jMenou. Perna &: Hernquist 



2001 



Eksi fc Alpar 



s used 



2003|) 



employ a gradual decrease of the mass-flow rate, M, at late ages, leading to a population 

00 s contrary to the c utoff period 



2002 ). Physically, 



of old observable sources with periods extending to 
P = 15 s inferred from an analysis of period statistics (iPsaltis fc Miller! 
as a disk evolves to decreasing M, it will cool starting from the outer disk. For a while 
irradiation from the neutron star may keep the disk sufficiently hot. As M decreases even 
irradiation will not sustain high temperatures. One expects that eventually the ionization 
fraction in disk matter will be so low that the magneto-rotational instability (MRI) (Balbus 
& Hawley 1991) which is believed to provide the turbulent disk viscosity should fail to 
operate. The disk in late phases should make a transition to an effectively neutral, passive 
state starting from the outer disk and leading to a cutoff in M and luminosity. We discuss 
the results of model calculations in 5 4 and summarize our conclusions in 8 5. 



2. EVOLUTION OF VISCOUS THIN DISKS IN TERMS OF INITIAL 

MASS AND ANGULAR MOMENTUM 



The evolution of the surface mass density X in a viscous thin disk (jPringle 
described by the diffusion equation 



1981) is 



( ILiist 



d£ _ 3 d 
dt r dr 



a/2 



— I 

dr 



19521 ) where v is the turbulent viscosity. For opacities of the for m k = Kop a T one can 



solve the disk structure equations to obtain viscosity in the form (see 

v = Cr p Z q 



Cannizzo et al. 



1990) 



where C, p and q are constants given by 



C 



Ql 6-26+a 



V32as B ; 



6-26+a 



6-2b+a 



+ 1 



(2) 



(3) 
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V 

(I 



-3 



a + 1 



1 



6-2b + a 2 
2 (a + 2) 



6 - 2b + a 



(4) 
(5) 



Here M* is the mass of the neutron star, a is the viscosity parameter ( jShakura fc Sunyaev 
19731 ). Jx is the mean molecular weight, k-Q is Boltzmann constant, m p is the mass of the 
proton, and (Jsb is the Stephan-Boltzmann constant. 

In general q ^ and v depends on the surface mass d ensity and the diffusion equation 



Pringld (119741 ) corresponding to 



([I]) is nonlinear. Three self-similar solutions were given by 
different boundary conditions. The solution of interest here is the one corresponding to a 
disk of constant total angular momentum and decreasing mass. This analytical model is 
a useful guide for the evolutionary calculation for the neutron star's spin, since the total 
angular momentum transfer from the neutron star to the disk is ne gligible c ompa red to the 



Pringld (119741 ) diverge 



angular momentum of the disk. The original solutions as given by 
at t = 0. As the diffusion equation is symmetric with respect to translations in time, we 
can shift the origin of time by t — > t + to for all occurrences of t in the solutions. With this 
modification the solution becomes 



E 



k 1 + 



t 

To 



5q-2p+4 



-Rout (t) 



V 
9 + 1 



1 



-Rout (t) 



2--2- 

9 + 1 



1/9 



where 



Rent (t) = r [ 1 + - 



(6) 



(7) 



is the outer boundary of the disk and 

k = 



(4q -2p + 4) (5q - 2p + 4) 

In these equations E , to an d ^o are constants introduced in the non-dimensionalization 
procedure and will be defined below. 
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Using the solution in equation ([6]), the total angular momentum of the disk is 
J d = J r 2 n K ll ■ 2ixrdr = 7 1 47rr 2 l £ V / GMr 



a constant. Here 



7i = k— — ^—mPv 



(9) 



(10) 



Aq-2p + A q Aq - 2p + 4 
is a constant where /3(k, I) = T(k)T(l)/T(k + /) is the beta function. The total mass of the 
disk 

i 

(11) 



E • 2-nrdr = M [ 1 + - 



M d = 

is a decreasing function where 

M = 47rr 2 S 72 (12) 
and r )i = kqj (Aq + 4 — 2p). Accordingly, the mass flow rate at the inner disk is 

(a - 1)M / t V" 

where a = 1 + l/(5g — 2p + 4). For a disk of given initial mass M and angular momentum 
Jo the other scale factors can be found through the following sequence: 



r 
So 

to 



(Jq/Mq) 2 72 , 

CM Si' 
M 



47rro7 2 

4r| 
3^o 



(14) 

(15) 
(16) 
(17) 



Cannizzo et al. 



(119901 ) showed that the self-similar solut ions of iPringld (11974J ) 



corresponds to the intermediate asymptotic ( IBarenblatt 



19961 ) behavior of a disk with 



the no- viscous-tor que boundary condition at the inner radius ensured by E(i2 in ,£) = 0. 
In Figure [T], we compare the analytical solution given in equation ([6]) with the numerical 
solution of the diffusion equation with a Gaussian initial distribution discretized on 2000 
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grids over a numerical domain extending from 10 8 cm to 10 13 cm. The values of constants 
are given in the Table [TJ For comparison, we have taken the initial mass and angular 
momentum of the disk to be the same in both solutions. We see that the numerical 
solution settles onto the analytical solution after the initial transient stage. The mismatch 
near the inner boundary, in the left panel, is a characteristic issue with the self-similar 
solutions which usually can not accommodate boundary conditions. Nevertheless, this is 
not important for calculating the mass flow rate at the inner radius. The panel on the right 
shows the mass flow rate at the inner disk which agrees well with the numerical solution 
after the brief transient stage. 



EDITOR: PLACE FIGURE Q] HERE. 



The early papers studying fallback disks for AXPs (IChatterjee. Hernquist fc Narayan 



2000 



Eksi fc Alpar 



2003 



Perna. Hernquist fc Narayan 



20001 ) employed the dynamical 



timescale (~ 1 ms on t he su rface of the neutron star) for t in equation ([TBI . Numerical 



calculations of 



Pringld (I199ll ) show that power-law evolution starts in a viscous time-scale. 



Comparing analytical and numerical solutions, we find that the characteristic viscous 
timescale to of the disk is determined by the initial mass and the total angular momentum 
of the disk. The equations (1171) can be combined to give 

4(47r 72 r( 7 2/7i) 4g - 2 p+4 Jo q ~ 2p+4 
3C(GM) 2 i-p+ 2 

For an electron scattering dominated disk, 

Jo 



Mi 



5g-2p+4 ' 



1* 







to ^ lOy 



lO^CH^S" 1 



7/3 



-2/3 



(19) 



1O" 4 M _ 

showing that the viscous timescale in the disk is very sensitive to the average value of 
specific angular momentum j in the disk. Indeed, to increases to ~ 2000 y if jo ~ 10 20 cm 2 
s _1 , a realistic range according to our calculations in § 3. 
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For a given initial mass of the disk, setting to corresponds to setting the initial angular 
momentum of the disk as the fourth parameter in addition to the initial mass of the disk 
together with the magnetic dipole moment and the initial period of the neutron star. 

The analytical solution presented in this section does not address transient behavior 
due to thermal- viscous instabilities and/or the variation of the active disk radius. In § 3, 
we will explain how these instabilities are introduced into the framework. 

We use the analytical model as a guide to the time consuming numerical calculations, 
for the choice of initial conditions and to exhibit the resulting neutron star luminosities as 
shown in Figure [2j This Figure shows the iso-luminosity lines of fallback disks at the age 
of 10 3 y and 10 4 y, respectively, depending on the initial mass and initial specific angular 
momentum of the disk. It is seen that there is a wide band of this parameter space that can 
lead to a disk evolution with L = 10~ 3 — 10~ 2 Le, the observed AXP luminosity range, at 
the age of 10 3 - 10 4 y. 



EDITOR: PLACE FIGURE [2] HERE. 



In the following section, for our numerical simulations, we will choose an initial disk 
density distribution corresponding to a specific angular momentum of ~ 9.7 x 10 19 cm 2 
s _1 which leads to AXP luminosities with an initial disk mass of 10~ 4 — 1O _5 M as can be 
seen from Figure [2J These values lead to viscous timescales of to ~ 2000 y and to ~ 9000 
y, respectively. These values for to are orders of magnitude greater than the dynamic al 



timescale employed in earlier models (lEksi fc Alpar 



2003 



Menou. Perna fc Hernquist 



200ll ). and lead to a slower decline of the disk mass and M. The inner radius of the disk, 
determined by M, will also be evolving on this long timescale in which the period of the 
neutron star evolves. This leads to an evolutionary picture very different from the earlier 
studies. In earlier work, with ad hoc choices of short t , the inner radius of the disk was 
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seen to move outwards very rapidly as the disk evolved. The inner disk radius could be 
kept inside the light cylinder only if the neutron star was subject to large torques, i.e. if it 
had large magnetic fields (0.5 — 1 x 10 13 G) so that the star spun-down sufficiently fast and 
the light cylinder expanded rapidly enough. The torque equilibrium for such high magnetic 
fields requires high mass flow rates. This led to high luminosities in model calculations. To 
avoid this problem, a rather efficient propeller was required allowing only ~ 10~ 2 of the mass 
inflow to accrete onto the neutron star. In the present model the realistic slow evolution of 
the disk allows neutron stars with smaller magnetic fields to become AXP/SGRs. 



We solve the dis 



MODEL CALCULATIONS 



£ diffusion equation (CO) for the surface density E of the disk (sec e.g. 



Frank. King fc Raine 



(120021 )) to examine the evolution of the X-ray luminosity and the 
spin period of the neutron star together. We take the initial surface density distribution as 
E oc r _A with A = 3/4 which is the expected surface- density profile of a thin, steady-state 
accretion disk. Independently from the actual initial mass distribution of a fallback disk, 
the disk's subsequent evolution quickly leads to an extended thin disk. 

Opacities depend on the details of the disk composition that are not well known for 
fallback disks. As in our earlier work on enhancement light curves of AXPs and SGRs, we 
use electron-scattering opacity in our calculations. Using bound-free opacities does not 
significantly affect our results. In any case the uncertainty in opacity is folded in with the 
uncertainty inherent in the a parameter of the viscosity. We solve the disk diffusion equatio n 



by em ploying the a prescription of the kinematic viscosity v = ac s h ( jShakura fc Sunyaev 



19731 ) where c s is the local sound speed and h is the pressure scale height of the disk. 

What should be the value of the a parameter to be used in our numerical calculations? 
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To answer this question, we pursue the results of our earlier work on the enhancement 
phases of persistent and transient AXP/SGRs. X-ray and infrared enhancement light 
curves of both persistent and transient AXP/SGRs can be explained by the evolution of 
th e disk after the inner disk matter has been swept outwards by a soft gamma-ray flare 



( lErtan fc Alpar 



2003 



Ertan. Gogiis fc Alpar 



2006 



Ertan fc Erkut 



20081 ). There is an 



important difference between the X-ray outburst light curves of persistent and transient 
AXPs and SGRs. For the persistent AXP/SGRs, model fits are obtained by using a single 
viscosity parameter (a ~ 0.1) that is similar to those employed in modeling the hot states 
of soft X-ray transients (SXTs) and dwarf novae (DN) in their outburst phases. This 
probably indicates that for persistent AXP/SGRs, the inner part of the fallback disk within 
a sufficiently large radius remains in the same viscosity state during the enhancement. 
Fitting the X-ray and IR enhancement light curves from persistent AXP/SGRs does 
not give information about whether the disk is in a colder (lower) viscosity state at the 
outer regions of a fallback disk. Help comes from the transient AXP/SGRs whose low 
luminosity quiescent phases and transient bright phases indicate transitions between 
different viscosity states of the active disk. The light curves of the transient sources can be 



reproduced if the disk u ndergoes a thermal-v i scous disk instabi 



T crit ~ 1000 - 2000 K (Ertan fc Erkut 



2008 



Cahskan &: Ertan 



ity at a critical temperature 



20091 ). As temperatures 



drop below T crit the viscosity parameter changes from ahot — 0.1 to a co id — 0.025 — 0.035 
in exactly the same way applied in the disk instability models of SXTs and DN for wh ich 



10 4 K corresponds to the ionization temperature of hydrogen (see e.g. 



Lasota 



(120011 ) for a review of the disk instability model). This difference between transient and 
persistent AXPs and SGRs could be attributed to their different accretion rates prior to 
the enhancement. In the long term, it is the cold state viscosity that drives the mass-flow 
rate from the outer to the inner disk in quiescence. The viscosity parameter cthot of the hot 
inner disk does not significantly change the long-term evolution of the disk. In all model 



13 



calculations, we take a co id = 0.03, «hot = 0.1, and T crit = 1000 K. We note that in both hot 
and cold states the disk is active, providing mass inflow with turbulent viscosities. The disk 
instabilities at T crit ~ 1000 K do not affect the long-term evolution. Using only a single 
Oi = a co id, without introducing a T crit in the model, the same results obtained for the long 
term evolution of the disk are almost the same as in the case of a two-state disk. 

Our main concern here is on how the long-term evolutions of AXPs and SGRs are 
affected if the disk becomes passive below some temperature T p which is much lower than 
and should not be confused with T crit . Termination of the disk's active lifetime by transition 
to a final passive phase at very low temperatures T p seems to characterize the M and P 
evolution of the AXPs and SGRs in the long term. 

The outer disk of a neutron star in X-ray binaries are cut by tidal forces. There are 
no such forces that can constrain the extension of a fallback disk around AXP/SGRs. 
Where is the outer radius r out of a fallback disk? By this r out we mean the radius beyond 
which active viscous mass transport ceases. A strong possibility is that the outermost 
disk could be passive at very low temperatures. To put the question in other words, what 
is the minimum temperature belo w which disk regions ent er a passive (non-turbulent) 



low viscosity regime? The work by llnutsuka fc Sand (120051 ) shows that the viscosities are 



)5) 



expected to be turbulent by the magneto-rotational instability (MRI) (Balbus Sz Hawley 
1991), even at temperatures ~ 300 K. In agreement with this result, in the present work, 
we show that a rather low minimum active disk temperature T p ~ 100 K, below which 
the viscosities are not turbulent, can explain the period and X-ray luminosity evolution 
of AXPs and SGRs consistently with the observed period, period-derivative and X-ray 
luminosity distributions of these sources. 

The minimum active disk temperature, T p , should not be confu sed with the cri t ical 



temperature, T crit , leading to a disk instability that was proposed by 



Ertan & Erkut 



(120081 ) 
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as the mechanism responsible for the X-ray enhancement of the transient AXP XTE 
J1810-197. Unlike T p , T crit has no significant effect on the long-term evolution of a fallback 
disk. It determines the transition between cold and hot state viscosities that governs the 
mass-flow rate, M, in the active parts of the disk. 

The position of r out is not important in simulating the short-term events like 
enhancements lasting from a few months to a few years as long as r out is taken sufficiently 
large in order to avoid numerical effects due to artificial cutoff of the disk. Numerically the 
outer disk must be cut somewhere. A constant r out is not realistic for a fallback disk and 
causes numerical effects on the model M evolution in the long term. The natural transition 
to a passive disk, occurring at T = T p , determines r out as the radius beyond which T < T p . 
Thus, r out is a dynamical quantity, moving from the outer disk inward as the disk ages and 
cools. 

From the model fits to the X-ray enhancement data of the transient AXPs and SGRs, 
we find that fallback disks are likely to be active even at temperatures about 1000 K. At 
much lower temperatures it is plausible to expect that a fallback disk gradually becomes 
passive starting from the outer radius which has the lowest temperature along the disk. 
At radii larger than some transition radius r p , at temperatures lower than a critical 
temperature T p , the outermost disk becomes passive. As the mass accretion rate decreases 
with time, leading to decreasing X-ray irradiation strength, the passive outermost disk 
grows (r p moves inwards). The transition temperature T p below which a fallback disk 
becomes locally passive is of crucial interest for the late evolution of the disk. 

We keep the inner disk radius r in constant. In the numerical model, separati o n of t he 



radial grid points decreases with decreasing radial distance r (see lErtan fc Alparl (120031 )). 
Narrower radial grid separation requires smaller viscous time step. Therefore, time steps 
for the calculations along the radial grids become smaller with decreasing inner disk radius. 
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Due to these numerical reasons, we take rj n = 10 10 cm which is larger than the expected 
inner disk radius of a fallback disk around an AXP with a conventional dipole magnetic 
field of strength ~ 10 12 G on the surface of the neutron star. This simplification does 
not affect the mass- flow rate at the inner disk, M- m . The realistic inner disk radius of the 
disk is expected to be close to the Alfven radius, r& = M in 2 ^ 7 (GM^ 1 ^ 7 fit^ 7 , where /x* is 
the dipole magnetic moment and M* is the mass of the neutron star. We calculate the 
disk torque acting on the neutron star through the interaction of the inner disk with the 



magnetic dipole field of the neutron star usi ng p,* = 10 30 G toget 



rer with the values of 



M in and ta calculated at each time step (see lErtan fc Erkutl (120081 ) for details). Thus, the 



fastness parameter and the torque are calculated dynamically, using the current r& value 
at each step. The fixed choice of r in , for numerical reasons, therefore does not influence the 
torque evolution either. 

The torque model is based on the angular momentum exchange between the stellar 
magnetosphere and the disk matter that penetrates through the effective magnetic coupling 
radius r\ in the fast rotator regime with lu* > 1, where lu* = (Va/V co ) 3 / 2 is the fastness 
parameter. We assume that the effective magnetic coupling radius is given by the Alfven 

1 /3 

radius ta- Here, the co-rotation radius, r co = (GM^/Vtl) , represents the critical radius 
where the disk matter rotates with the same speed as the neutron star. This is true if the 
rotation of the disk is purely Keplerian. In reality, the actual rotation rate of the disk 
matter slightly exceeds the local Keplerian value at each radius r between the effective 
co-rotation radius and the Alfven radius. This is the result of the transfer of angular 
momentum from the neutron star to the inner portion of the disk threaded by the stellar 
magnetic field lines over the range r co < r < In this regime of spin-down with accretion, 
the magnetic stresses gradually force the inner disk matter to co-rotate with the neutron 
star within a transition region in the disk where the rotation rate Q changes from nearly 
Keplerian values in the outer parts to co-rotation with the magnetosphere at the inner 
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boundary of the disk. 



For a magnetically coupled disk-neutron star system, the integration of magnetic 
stresses over the magnetically th readed region o f the d isk yields the effective magnetospheric 



Erkut fc Afparl (120041 ) ) . In the mildly fast rotator regime 



torque on the neutron star (see 
(cu* > 1), assuming most of the incoming disk matter, instead of being propelled out of 
the system, accretes onto the neutron star, the effective magnetospheric torque can be 
estimated as 



N* ~ - / r 2 B+B z dr = yMi/ GMta (l - 

J Tret 



(20) 



where B^ = 7^i? 2 is the azimuthal component of the magnetic field on the surface of 
the disk, 7^ > is the azimuthal pitch of order unity for r > r co , and B z ~ —fi^/r 3 is 



the vertical c omponent of the stel 



r C o < r < ta (jErtan fc Erkut 



200Sh . 



ar dipole field threading the inner disk over the range 



Ertan&Erkut 



(J2008|) to 



The disk torque in equation (|20|) has been recently used by 
explain the spin-frequency evolution of the transient AXP XTE J1810 — 197 during the 
long-term post-burst fading of its X-ray flux. The disk model has also acco unted for the 



X-ray light curve of the AXP XTE J1810 — 197 during the same period (see 



Ertan&Erkut 



( 120081 )). The model fits to both the X-ray light curve and spin- frequency data are consistent 
with the regime in which the neutron star spins down while accreting most of the matter 
flowing in through the disk. When the inner-disk radius is outside the light cylinder, we 
assume in the torque calculations that the inner disk extends down to the light cylinder and 
that all the mass is propelled out of the system. In the model calculations, this happens 
during the early evolutionary phases, and these assumptions do not change the subsequent 
evolution curves of the model sources. 

As in the case of the AXP XTE J1810 — 197, X-ray luminosity changes in the transient 
AXPs and SGRs cover a broad range from ~ 10 35 erg s" 1 down to ~ 10 33 erg s _1 . Such 
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a wide range of X-ray luminosities, in a given transient source, encompasses both the 
relatively high accretion regime (~ 10 14 — 10 16 g s _1 ) of the persistent AXPs and SGRs and 
the low a ccretion regime ( S, 10 1 3 g s _1 ) in the quiescent phase of their transient cousins. As 



shown by 



Ertan fc Erkut 



(120081 ). the disk torque in equation (120]) describes well the spin 
evolution of the neutron star during the long-term X-ray enhancement of the transient AXP 
XTE J1810 — 197 while the X-ray flux changes by up to two orders of magnitude. The 
dependence of the torque on the fastness parameter, that is, iV* oc 1 — u%, accounts for the 
spin-frequency evolution of the relatively fast rotating (a?* > 1) neutron stars for which the 
long-term decay of the X-ray lu minosities can be expla ined by a decrease in the mass-inflow 



rate in the inner disk (see, e.g., 



Ertan fc Erkut 



(120081 )). The same torque, on the other 



hand, behaves as A* oc 1 — cj* for the n eutron stars that are close to rotat ional equilibrium 



Ertan. Erkut. Eksi & Alparl (12007m . Whether they 



with the disk, i.e., for uo* ~ 1 (see, e.g., 
are close to rotational equilibrium or not, the torque in equation (120]) is applicable to both 
the persistent and transient sources and thus will be employed in our calculations for the 
evolution of the spin-period and X-ray luminosities of the AXPs and SGRs. 



4. RESULTS AND DISCUSSION 

In our calculations, we see that for a given initial disk mass, M , and a minimum active 
disk temperature T p , there is a minimum initial period Po,min such that sources with initial 
period Pq < Po,min cannot evolve into the AXP phase. For the sources with P > Po,min, 
time for the onset of accretion depends on Pq. The neutron stars having longer Pq start 
accreting matter earlier. For a large range of M , all the model sources that can enter the 
regime of accretion with spin-down converge into the AXP/SGR phase. When the initial 
mass of the disk is too large, above a certain M , X-ray luminosities remain higher than 10 36 
erg s -1 (or equivalently M > 10 16 g s^ 1 ) for sufficiently long time intervals, which would 
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make sources with L x > 10 36 erg s 1 . Since, such luminuous sources are not observed, there 
is probably an upper limit to M values. 

In the first set of simulations (Figs. EH6]), we set T p = 100 K. For different values 
of initial total disk mass Mo, we obtain M and period evolution curves. In the figures 
we have indicated the lines for M = 10 14 g s _1 and M = 10 16 g s _1 to delineate the 
M interval corresponding to the observed AXP/SGRs. In Figure [3], it is seen that for 
M = 1 x 1O" 5 M only the sources with Pq > 70 s can evolve to the AXP phase, and those 
with lower P are likely to become, and remain as radio pulsars throughout their lifetimes. 
Model sources with initial periods 80, 100, and 1000 ms start accretion at different times, 
while all of them evolve into the AXP phase. The minimum initial period required for the 
source to eventually enter the AXP phase, -Po.mm, is correlated with the initial mass of the 
fallback disk; Po,min decreases with increasing M . For M values of 4 x 10~ 5 , 1.0 x 10~ 4 
and 3.5 x 10 _4 M Q we obtain minimum initial periods around 30 , 15 and less than 10 ms 
respectively (Figs. |4]-[6]). For a large range of Mo values, Figures [3] - [6] clearly show that 
the sources with P > Po,min converge to the AXP/SGR phase and more significantly, their 
period evolution converges to, and remains in the period range of AXP/SGRs as their 
luminosities decrease below minimum AXP luminosities. Given the observed range of X-ray 
luminosities, the model curves seen in Figure 6 are likely to represent the evolution with an 
extremely high M that is probably either unrealistic or not common. It is remarkable that, 
even in this case, M and P curves tend to meet in the AXP/SGR phase. 

How do these results change if we set T p = 0, that is, if the disk remains active for 



all temperatures? This is actually an impli cit assumption in all earlier ana- 



using self-similar evolution of fallback disks (IMenou. Perna fc Hernquist 



2001 



ytical work 



Eksi & Alpar 



20031 ). For this test, we set M = 1 x 10 4 M Q (same as for the models given in Fig. [5]) and 



T p = without changing the other parameters. We repeat the calculations for different 
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initial periods. The resultant model curves are presented in Figure [7J It is seen in these 
figures that model sources starting with different neutron star initial periods do settle into 
the AXP/SGR period and luminosity range. As expected, there is a crucial difference 
between the models with T p = K and T p = 100 K in that the cutoffs in M curves for 
T p = 100 K are not seen in M evolutions with T p = 0. This means that the mass-inflow 
rate remains large, above or slightly below 10 16 g s -1 , for 10 5 years or more. Thus, without 
a cutoff in disk activity, agreement with the observed AXP/SGR luminosity range would 
require an explanation. Ek§i & Alpar (2003) suggested that a large fraction of mass flowing 
in through the disk is not actually accreted onto the neutron star, requiring rather efficient 
propeller activity. 

We have calculated more models for a better comparison of evolution with and without 
the disk becoming eventually passive. First, we take P — 300 ms, T p = and repeat the 
calculations for different M values (Fig. [8]). (We note that the model sources with low 
Mq values (curves 5 and 6 in Figure [8]) that are seen to evolve towards the AXP phase 
with Pq = 300 ms, would become radio pulsars for initial periods Pq < 100 s). This shows 
that sources with a large range of initial periods evolve into the AXP range once they have 
started to accrete matter. The model curves in Figure [8] correspond to different initial disk 
masses ranging from 7.5 x 10 -7 to 1 x 10~ 4 M Q . All these illustrative models remain within 
the range of AXP/SGR luminosities for more than 10 5 y. During most of this epoch, the 
sources lie in the AXP/SGR period range as well. Nevertheless, if this model were a correct 
representation of these systems we would expect to see many more, mostly older, AXPs and 
SGRs with similar luminosities and with periods larger than the observed range of period 
clustering. This is a serious difficulty for fallback disk models without an eventual cutoff for 
M. 

Next, we do the same exercise with the same Pq and the same set of M Q values, but for 
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T p = 100 K (see Figure [9]). The model sources that enter the period range of AXP/SGRs 
follow almost the same period-evolution curve. They enter the AXP phase and remain 
there for a long time interval lasting from a few xlO 3 y to a few xlO 4 y depending on the 
initial mass of the fallback disk. Subsequently, periods of the sources tend to increase with 
about the same rate to beyond the maximum observed AXP periods. By comparing panels 
a and b in Figure M, we see that the sudden turn down in X-ray luminosity renders the 
sources unobservable at an age of about a few 10 4 y. This terminates the observable period 
evolution in an epoch when the periods are clustered in the 2 — 12 s range, in agreement 
with observations. 

For given initial disk mass Mq, how does the evolution depend on T p ? To illustrate, 
we compare the M and period evolution with T p values of 0, 60, 80, 100 and 200 K for 
Mo = 1 x 10 4 M Q and Pq = 30 ms. For this Mq, the model sources having T p greater 
than about 250 K cannot evolve into the AXP phase. This is because higher T p means 
earlier decrease in mass-flow rate from the outer to the inner disk, and thus the inner disk 
cannot penetrate into the light cylinder (see Figure [T0i) . Taking also the earlier model fits 
with different disk masses into account, the results obtained with T p ~ 100 K are in good 
agreement with observations of AXPs and SGRs. 

In our model calculations, we take B = 10 12 G. For larger dipole magnetic fields 
of the neutron star the eventual, approximately equilibrium periods reached in the late 
evolutionary epochs would be longer in proportion to B 6 ^ 7 . To maintain agreement with the 
observed period clustering requires earlier luminosity cutoff. In the model, for larger dipole 
magnetic fields of the neutron star up to about 1 x 10 13 G, it is still possible to reproduce 
AXPs by increasing T p and Mq values accordingly. Increasing B results in a narrower initial 
Mq range that can give AXPs. This is because increasing field strength requires higher 
Mq for the inner disk to be able to penetrate the light-cylinder. For the field strengths 
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above ~ 10 13 G, model sources cannot enter the AXP phase, both the luminosity and the 
period derivative remain well above the observed values while the neutron star is tracing 
the AXP period range. Note that the .B-field discussed here is the neutron star surface 
value of the dipole component, as this is the long range component of the magnetic field 
that determines the torques and interactions with the disk. The magnetar fields of AXPs 
and SGRs are likely to be in higher multipoles in view of the crustal field amplification 
mechanisms responsible for magnetar strength fields. It is the dipole fields that must ha ve 



magnitudes less than about 10 13 G. Recent work on AXP surface spectra (iGiiver et al. 



2007 



20081 ) has produced estimates of the total surface magnetic field strengths, in AXPs 
XTE 1810-197 and 4U 0142+61. The total surface fields inferred were ~3x 10 14 G and 
~ 5 x 10 14 G respectively, larger than the B < 10 13 G requirement of our model for the 
dipole component alone. Our result also implies that some of the young neutron stars that 
have high dipole fields B > 10 13 G are not likely to evolve into the AXP phase even if there 
are fallback disks around them. 

This paper does not address the relative strengths of dipole and higher multipole 
components on the basis of neutron star physics. We do think, however, that dipole fields 
weaker than the higher multipole fields are plausible for magnetar models. The large fields 
in the star's crust are amplified, rearranged in crust breaking events and dissipated in the 
neutron star crust. Magnetic stresses are here in interaction with crustal crystal stresses, 
which have critical strains, dislocation and domain wall distributions that introduce local 
scales to the problem. In the neutron star crust with very efficient screening, the crystal is 
close to a Coulomb lattice. So the relevant scales of microcrystalline alignment can be much 
longer than those in terrestrial crystals, but still they are small compared to the global 
size of the neutron star. Field generation and reconnection in co upling with discr ete crust 



breaking effects will plausibly put energy into higher multipoles ( IRuderman 



199J ). 
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In P, P and timing signatures, there is a continuity between AXPs, SGRs and the 
radio pulsars. If the torque causing the spin-down is the dipole spin-down torque for an 
isolated neutron star, than the dipole magnetic field is a function of P and P ( assuming 



all neutron stars in the s ample have standard values o 



of this were discussed by 



Rt/ I\ possible effects of variations 



Guseinov. Ankay &: Tagieval (120051 ) ) . Then the neutron stars with 



similar positions in the P — P diagram would be continuous in -B^dipoie also, leaving a puzzle 
as to why some are AXPs and SGRs and others not. Furthermore, if the dipole spindown 
torque is operating on all AXPs, SGRs and high -B* i di P oie pulsars, evolutionary connection 
between AXP/SGRs and high -B*,di P oie radio pulsars would require ongoing dipole field 
generation, as constant field evolutionary tracks for isolated neutron stars have slope -1 in 
the log P — log P diagram. 

Presence or absence of fallback disks and their properties bring in a new parameter, 
not determined by P and P alone, into play. The dipole fields inferred by using fallback 
disk model spindown torques are somewhat lower than the dipole fields inferred using the 
isolated rotating dipole model. The transformation from P and P to -B*,dipoie involves 
the mass inflow rate M and fastness parameter, whose values may be constrained by 
observations of the X-ray luminosity for some sources. Discussion of the P — P diagram in 
the light of this transformation will be the subject of our subsequent work. In any case, 
sources located close to each other in the P — P diagram are not likely to have similar 
values of -B*,dipoie if their spindown is due to fallback disk torques. Inferred dipole fields of a 
radio pulsar and an AXP lying close to each other in the P — P diagram will be different 
also in the case that the radio pulsar does not have a fallback disk, so it has a relatively 
higher -B^dipoie inferred by the dipole spindown torque, compared to the AXP which will 
have a relatively lower -B^dipoie inferred by the disk spindown torque. 

If there are radio pulsars with fallback disks, their surface dipole magnetic fields, 
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calculated correctly with the disk torques are likely to be in the -B*,di P oie ~ 10 12 G range. 
These can evolve to become AXPs. Our results show that pulsars with -B* idipo i e > 10 13 G 
are not likely to become AXPs. 

Why do the high PP radio pulsars (except one case) not show bursts (do not have 
higher multipoles? After the dipole fields are calculated with the appropriate disk torque 
for those neutron stars with fallback disks, there may remain still some radio pulsars with 
surface dipole fields comparable to those of AXPs. There is also t he one example (at 
present) of a bursting radio pulsar PSR J1846-0258 in SNR Kes 57 (IGavriil et al.lhoosh . So 
the question arises as to what makes an exceptional bursting radio pulsar. If as we claim 
the burst behavior requires magnetar-strength fields in higher multipoles due to a history of 
crust breaking, than the usual, non bursting radio pulsars (in general) do not have strong 
higher multipoles even if some of them may have rather strong surface dipole fields. A likely 
explanation is that magnetar field evolution starts with the initial total magnetic moment 
of the neutron star, and is driven by the initial field in the neutron star core. Suppose 
a magnetar field in higher multipoles would obtain if the initial total field in the neutron 
star's core was higher than some critical value, enough to affect crust breaking (say by the 
magnetic flux lines extending through the superconducting core). The high dipole field 
radio pulsars started with fields near or just below the threshold for generating higher 
multipole magnetar fields: thus the dipole field in these pulsars, like in other radio pulsars 
is roughly the total field. The one exception, PSR J 1846-0258, may be the one marginally 
above threshold in its initial neutron star core magnetic moment. 

Finally, we compare the P values of the model sources during their AXP/SGR phase 
with observations. For all models with different Mq and Pq reaching the AXP/SGR phase, 
we see that P values remain in the range from a few xl0~ 13 to ~ 10~ n s s" 1 . This 
is consistent with the observed period derivatives of these sources. We present P and 
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fastness variations for the same illustrative model seen in Figure [TU1 with T p = and 100 K 
(Fig. ITTk). The corresponding evolution of the fastness parameter, a;*, for the same model 
is given in Figure [TTb. The maximum of a;* at t ~ 2000 y corresponds to the time at which 
the inner disk penetrates the light cylinder. Before t ~ 2000 y, the inner disk remains 
outside the light cylinder. In this early phase of evolution, during which we take the inner 
disk radius equal to the light cylinder radius, the strength of the disk torque increases with 
increasing light-cylinder radius caused by rapid spin-down of the neutron star. For this 
example, after t ~ 2000 y, the inner disk remains inside the light- cylinder and the model 
sources evolve first to the so called tracking phase along which uj* stays close to unity, that 
is, inner disk comes close to the co-rotation radius. In Figure [TTb . it is seen that cu* of the 
source with T p = 100 K turns up from the tracking phase at about a few 10 4 y and starts to 
increase, while the model source with T p = remains in the tracking phase. This is because 
the mass-flow rate for T p = 100 K decreases faster than that for T p = K. The resultant 
rapid increase in the inner disk radius for T p = 100 K causes a;* to increase to higher values. 
The source is likely to become a strong propeller. A fraction of the mass would be expelled 
from the system, instead of accreting onto the neutron star. Qualitatively this would further 
decrease the observed luminosity and mak e the luminosity turn -off even sharper. Work on 



the luminosity changes of transient AXPs (lErtan Erkut 



20081 ) shows that mass expulsion 



does not yield a se lf-consistent mod e 



for 



uminosities down to ~ 10 33 erg s 1 , in conclusion 



(120041 ) . The effect needs to be investigated at lower M 



in agreement with iRappaport et al. 
(and higher u*). 

There are several possibilities for the radio emission in AXPs with fallback disks. 
One of them is the conventional radio-pulsar emission mechanism. It is possible that this 
mechanism could work in the presence of accretion, with the strong magnetic field confining 
the a ccretion column. This could shiel d and allow the radio emission along the open field 



lines (jTrumper. Ertan. fc Kylafis 



2009L 



The accretion column is likely to have a thin 



-25 - 



cylindrical-shell geometry enclosed by the closed field lines with a radius that depends on 
the inner disk radius, where the inner disk matter couples to the closed field lines. Accretion 
energy is injected into the neutron star through the base of the column and radiated from 
the entire surface of the neutron star in the soft X-ray band, while the column could act as 
the site of the beamed high-energy emission. In neutron stars with long spin-periods, like 
the AXPs, the radio emission is expected to be strongly beamed with a full beaming angle 
of several degrees along the open field lines. This indicates that the radio emission efficiency 
in these systems is much lower than that estimated by assuming isotropic radio emission. 
Small beaming angles could explain why the pulsed radio emission is observed from only 2 
of the AXPs. 



A second alternative could be t 



re radio emission through 



disk-star radio emission mechanism (jCheng fc Ruderman 



1991 



gaps provided 



Ertan fc Cheng 



py th e 



200J). In 



this picture, the accretion torque, rather than the magnetic dipole torque, powers the gap. 

Another possibility for the radio emission could be related to charge acceleration by 
magnetar multi-pole fields close to the neutron star as suggested by Thompson (2008) if the 
site of the radio emission can be shielded from the accretion column. 



Comparing the various radio emission models with the radio and X-ray behavior of the 
transient AXPs XTE J1810-197 and IE 1547-54 will be the subject of our future work 



5. CONCLUSIONS 

By means of a series of numerical simulations of long-term evolution of the neutron stars 
interacting with fallback disks we have shown that: (1) The torque model iV* oc (1— uf) leads 
to the observed period clustering for a wide range of initial conditions provided the neutron 
star gets into and remains in an epoch of accretion with spin-down, which holds when the 
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inner edge of the accretion disk penetrates the light cylinder. (2) The correct identification 
of t yields long initial disk evolutionary timescale t ~ a few 10 3 y, without requiring 
excessive values of M and the disk mass. A wide range of initial conditions do lead to the 
epoch of accretion with spin-down. (3) The introduction of T p , the physically expected 
criterion for the disk to become neutral and passive starting from its outer regions, leads to 
evolutionary models remaining in the AXP/SGR X-ray luminosity range at the right ages 
when the period evolution has converged into the observed range of period clustering. We 
find the minimum temperature T p for an active disk to be T p ~ 100 K. This value is in 
agreement with the results of Inutsuka and Sano (2005), who concluded that turbulence 
works down to temperatures as low as ~300 K. The low values of T p indicates that the inner 
fallback disks are active up to a large radius greater than about 10 12 cm. (4) The period 
derivatives are also concur with observed values. (5) For dipole field strength > 10 13 G, we 
see that the model sources cannot evolve into the AXP phase consistently with observed 
luminosity, period and period-derivatives. Based on this result, we expect that sources 
with dipole fields of <; 10 13 G are not likely to become AXPs, even if they are born with 
fallback disks. (6) The transient AXPs and SGRs with X-ray luminosities L x ~ 10 33 erg 
s _1 in quiescence are likely to be older than the persistent sources. This is because fallback 
disks that have accretion rates corresponding to these luminosities are unlikely to penetrate 
the light cylinders of young neutron stars with small initial spin periods. The number 
of observable transient AXP/SGRs are expected to be comparable to those of persistent 
sources, as the models indicate that the transient AXP/SGRs are in the cutoff phase of the 
X-ray luminosities for a wide range of initial conditions (see Figure [9]) . 

These results could provide an opportunity to study the evolutionary connection 
between AXP/SGRs and other young neutron star popula t ions t hat could also be evolving 



with active fallback disks, as proposed originally by 



Alparl (1200 ll ). This will be the subject 



of our future work, in particular on comparisons with the central compact objects (CCOs) 
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and the dim isolated neutron stars (DINs). 
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Table 1. Values of some constants in two different opacity regimes 



Opacity 


c 


P <? 


a 


71 


72 


^cs ^ > ^bf 
/^bf --^ ^es 


\ GMa J \ /amp / 

4(^) 7 fe)"(GM)-A 


1 2/3 
15/14 3/7 


19/16 
5/4 


3.6851 x 1(T 4 
1.7030 x 10" 5 


6.2626 x 10~ 4 
3.3558 x 10~ 5 



Note. — The values of the constants in two different opacity regimes. Here, K es is the electron scattering 
opacity and n b j is the bound-free opacity. a ss is the Shakura Sunyaev viscosity parameter, fi is the mean 
molecular weight, M is the mass of the neutron star which we assume to be 1.4Mq. 



- 32 - 





-^Satoo-i • — - / ; 


Analytical - 

Numerical 

=0 years 
















"""^^i|S£s / 




































Fig. 1. — Comparison of numerical solution of the diffusion equation with the analytical 
solution given in equation [6j The initial surface mass density distribution is a gaussian 
E(r, 0) = £ max exp( — (r — r max ) 2 /(X 2 ) where r max = 2 x 10 10 cm, er = 4 x 10 9 cm and 
S max = 2.23 x 10 7 g cm -2 is chosen such that the initial mass of the disk is M = 1O _5 M . 
The corresponding initial angular momentum of the disk is Jo = 1-03 x 10 47 g cm 2 s _1 . The 
initial mass and angular momentum of the disk yield the following non-dimensionalization 
scales: S = 7.35 x 10 8 cm~ 2 , r = 5.86 x 10 10 cm and t = 5.74 x 10 2 y. Left Panel: 
Evolution of the surface mass density S. Right Panel: Evolution of the mass flow rate 
at the inner radius. 
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Fig. 2. — The X-ray luminosity of accretion from a fallback disk for a range of initial disk 
masses and initial specific angular momentum at t — 10 3 y (left panel) and at t — 10 4 y 
(right panel). 



-34- 




Fig. 3. — (Left panel) M evolution with initial mass M = 1 x 1O _5 M of the extended disk 
and with the transition temperature between active and passive viscosity states T p = 100 K. 
The models with the shown initial periods (in ms) start on the accretion phase of the same 
evolutionary track at the times indicated with dashed vertical lines. For this M , sources 
with P less than about ~ 70 ms are likely to evolve as radio pulsars (see the text for details). 
The two small spikes result from local viscous instabilities in the disk and are not numeric 
artifacts. (Right panel) Period evolution with the same parameters given in the left panel. 
Model curves with P values of 80, 100 and 1000 ms (noted on the left panel) are seen to 
enter the AXP/SGR phase. The sources with initial periods 40, 50 and 60 ms do not reach 
the spin-down with accretion regime. 



- 35 - 







50 


40 


\ : 








30 ■ ■ \ : 



10000 
Time (Years 



1 1 ' ' ; 














50 .... 


40y' 


30/' 
""' 15 










, , , , 10 





1000 

Time (Years 



Fig. 4. — (Left panel) Same model as given in Fig. [3] with Mo = 4 x 10~ 5 M Q , and initial 
periods 30, 40 and 50 ms. (Right panel) Period evolution for the same model with M curves 
seen in the left panel. With this initial disk mass, the model sources having initial periods 
10 and 15 ms cannot become AXP/SGR. 
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Fig. 5. — (Left panel) The same model as given in Fig. [3X with M = 1 x 10~ 4 M Q . Initial 
periods are 15, 20, 30 and 40 ms.; (Right panel) Period evolution corresponding to the M 
curves given in the left panel. The source having initial period of 10 ms cannot become an 
AXP for this initial mass, while all the other sources, comparing with left panel, are seen to 
evolve into AXP phase. 
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Fig. 6. — (Left panel) The same model as given in Fig. [3] except Mo = 3.5 x 1O~ 4 M . 
Initial periods are 10, 20 and 30 ms. This initial mass seems to be unlikely for fallback disks 
unless a large fraction of the inflowing disk matter is propelled out of the system for an 
initial evolutionary phase of a few xlO 4 y. Otherwise, we would expect to see sources with 
luminosities above the observed AXP luminosities and with periods around or less than the 
AXP periods (see the right panel). (Right panel) Period evolution corresponding to the M 
curves given in the left panel. 
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Fig. 7. — (Left panel) The only difference from the model given in Fig. [5] is T p = 0, that 
is, the entire disk is taken to be active throughout the calculations. The model M curves 
shown in the figure are obtained for different initial spin periods: 10, 15, 20, 30 and 40 ms. 
(Right panel) Period evolution for the same model model given in the left panel. 
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Fig. 8.— (Left panel) M evolution for T p = with different M values: 1 x 1(T 4 M (1), 
5.8 x 1O- 5 M (2), 2.5 x 1O~ 5 M (3), 3.5 x 1O" 6 M (4), 1.6 x 1CT 6 M (5), and 7.5 x 1(T 7 M 
(6). For a comparison of the model curves we take Pq = 300 ms for all the models (see 
the text). (Right panel) Period evolutions accompanying the M evolutions given in the 
left panel. Period clustering is also expected for T p = 0. There is an important difference 
between the M evolution curves of the models with T p = and T p = 100 K (see Figure [9]). 
For T p = 0, we would expect to see many AXPs and SGRs with ages greater than 10 5 y and 
with periods even larger than observed AXP and SGR periods (see the text for discussion). 
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Fig. 9. — (Left panel) The same models as given in Figure [H with T p = 100 K here. It 
is clearly seen that there is a natural cutoff for M and thus for X-ray luminosity. We see 
by comparing to right panel that statistical expectation of this model seems to be in good 
agreement with observations. (Right panel) Period evolution for the model curves given in 
the left panel. Model (6) with M = 7.5 x 1O _7 M is not present in the left panel, since its 
inner disk radius cannot penetrate into the light cylinder and the source evolves as a radio 
pulsar. It is striking that all the other sources which enter the spin-down with accretion 
phase also remain in a long lasting AXP/SGR phase. Furthermore, the model predicts that 
there will be an upper limit to the observed periods of AXP/SGRs in their present luminosity 
range because of the luminosity cutoff at the end of the evolution. 
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Fig. 10. — (Left panel) The same model as presented in Figure El but with different T p 
values which are seen on the curves in Kelvins. The initial period is 30 ms for all the model 
sources. With this initial period, the model sources having T p greater than about 250 K 
cannot become an AXP, since higher T p values lead to lower mass-flow rate from the outer 
to the inner disk. (Right panel) Period evolution of the same model sources given on the left 
panel. 




Fig. 11. — (Left panel) The period-derivative evolution with T p = and T p = 100 K. For 
both models, M = 1 x 10 _4 M Q and P = 30 ms. (Right panel) Evolution of the fastness 
parameter for the same model as given in the left panel 



